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MHD TIDAL WAVES ON A SPINNING MAGNETIC COMPACT STAR 
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ABSTRACT 

In an X-ray binary system, the companion star feeds the compact neutron star with plasma materials 
via accretions. The spinning neutron star is likely covered with a thin "magnetized ocean" and may 
support magnetohydrodynamic (MHD) tidal waves. While modulating the thermal properties of the 
ocean, MHD tidal waves periodically shake the base of the stellar magnetosphere that traps energetic 
particles, including radiating relativistic electrons. For a radio pulsar, MHD tidal waves in the stellar 
surface layer may modulate radio emission processes and leave indelible signatures on timescales different 
from the spin period. Accretion activities are capable of exciting these waves but may also obstruct or 
obscure their detections meanwhile. Under fortuitous conditions, MHD tidal waves might be detectable 
and offer valuable means to probe properties of the underlying neutron star. Similar situations may 
also occur for a cataclysmic variable - an accretion binary system that contains a rotating magnetic 
white dwarf. This Letter presents the theory for MHD tidal waves in the magnetized ocean of a rotating 
degenerate star and emphasizes their potential diagnostics in X-ray and radio emissions. 

Subject headings: MHD - stars: neutron - stars: oscillations - X-rays: stars - stars: binaries — stars: 
novae, cataclysmic variables 



1. INTRODUCTION 

For a magnetized spinning degenerate compact star cov- 
ered with a thin plasma layer referred to as an ocean, var- 
ious magnetohydrodynamic (MHD) tidal waves, excited 
by relevant processes, may exist. This scenario is plau- 
sible for an accretion binary system in which the Roche- 
overflow from a companion star accretes onto a spinning 
magnetized compact star. This Letter describes the theory 
of MHD tidal waves that propagate in a thin magnetized 
ocean of a spinning compact star. When combined with 
pertinent dynamic processes and emission mechanisms, 
such MHD tidal waves may leave indelible signatures of 
electromagnetic emissions in a wide spectrum ranging from 
the radio to GeV gamma rays and could provide key diag- 
nostics that might offer clues and hints of otherwise com- 
plex phenomena. In general, MHD tidal waves may be ex- 
cited gravitationally, thermally, electromagnetically, or by 
impacts of accreting materials. While accretion activities 
may stir up and sustain these waves, wave signals might be 
undesirably entangled with or overwhelmed by accretion 
noises. Once observationally established, these waves may 
be utilized to probe the structures in the surface layers of 
neutron stars and white dwarfs by extending or adapting 
the helioseismological techniques (Papaloizou & Pringlc 
1978; van Horn 1980; Saio 1982; Lou 2000, 2001b). Fur- 
thermore, damping properties of these waves might pro- 
vide information of the stellar constituent materials and 
diffusion processes (Lou 1995). 

2. THEORETICAL MODEL ANALYSIS AND RESULTS 

MHD tidal waves may exist on spinning magnetic de- 
generate stars in either isolation or binaries (Chanmugam 
1992; Phinney & Kulkarni 1994; Bildsten & Strohmayer 
1999), although they might have a better chance to be de- 



tected in high-mass and low-mass x-ray binaries (HMXBs 
and LMXBs) or cataclysmic variables. By strong grav- 
ity, the magnetized plasma ocean on a degenerate star of 
radius R is thin. As the ratio eR 3 flB/(GMm p c) of elec- 
tric force to surface gravity can be large for either neu- 
tron stars or white dwarfs in rotation (Goldreich & Julian 
1969), where e is the electron charge, f2 is the stellar angu- 
lar spin rate, B is the poloidal magnetic field, c is the speed 
of light, G is the gravitational constant, M is the stellar 
mass and m p is the proton mass, excess surface charges will 
be pulled upwards against gravity to form stellar magne- 
tospheres. It is presumed that the ocean and an extended 
magnetosphere coexist on a compact star, that the magne- 
tosphere maintains a pressure at the interface between the 
ocean and the magnetosphere much as the Earth's atmo- 
sphere does on the ocean (Gill 1982), and that magnetic 
fields (idealized as a split monople) are in and out radially 
across the thin ocean of thickness h. For a zonal magnetic 
field, the Alfvcn-Rossby waves have been studied earlier in 
the solar/stellar context (Gilman 1969; Lou 1987, 2000). 

For h/L <C 1 and a small Rossby number 1Z = U/(2QL) 
where U and L are typical horizontal velocity and spatial 
scales, the shallow-magnetofluid approach is valid. With a 
quasi-static vertical balance, the MHD tidal equations in 
spherical polar coordinates (r, 0, <fi) are: 
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where variables v$, v$, n and b r are the perturbations of 

^—component velocity, 0— component velocity, ocean sur- 
face elevation, and f— component of magnetic field, respec- 
tively, / = 2flcos9, g = GM/R 2 , B and p are the Coriolis 
parameter at colatitude 9, the surface gravity, the radial 
magnetic field, and the mean mass density, respectively, 
C A = B/iAirp) 1 / 2 is the Alfven speed, 6 = gr]+Bb r /(4wp) 
is proportional to the total pressure perturbation, and the 
displacement current effect is retained in case of Ca > c 
(Lou 1995). The energy conservation follows from eqns. 
(l)-(4) with the energy flux density v±[pgrj + Bb r /(4ir)] 
where v± = (vg,v<f,). By the exp[i(ujt + m<j>)] dependence 
with angular frequency u> and integer m, the resulting 
equation of Q(p) becomes 
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(gh + C 2 A )(l + C 2 A /c 2 ) (5) 

where p = cos0 and s = u(l + C\/c 2 )/(2Q). Variables 
vg, v<t„ n and b r may be explicitly expressed in terms of 
O(p). The 0(/i) solutions, regular at the two poles, in- 
volve the Hough functions (Hough 1897, 1898; Flattery 
1967; Longuet-Higgins 1968; Chapman & Lindzcn 1970). 
Given relevant physical parameters, all properties of MHD 
tidal waves can be derived computationally. 

Without magnetic field, eqn (5) is equivalent to eqns (7) 
and (8) of Bildsten et al. (1996) with their two parame- 
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q 2 = s~ 2 and A = R 2 uj 2 /(gh). In their "traditional 
approximation" , the thin compressible layer is disturbed 
adiabatically, the radial component of the Coriolis force 
is ignored, the radial displacement is much less than the 
horizontal displacement, and u> is much lower than the 
Brunt- Vaisala frequency. They focus on internal grav- 
ity modes of a rotating ocean. While internal magneto- 
gravity modes (Lou 1996) are excluded by invoking the 
quasi-static vertical balance, our incompressible formula- 
tion contains various MHD tidal waves in different regimes 
associated with the distinguished names of Kelvin (1879), 
Alfven (1942), Rossby (1939), and Poincare (1910), lead- 
ing to the acronym of KARP waves. 

For analytical simplicity, we hereafter invoke the 
(3— plane approximation (Rossby ct al. 1939) that re- 
tains the essential physics in the spherical geometry. In 
the local (3— plane approximation of eqns (1) — (4) with 
(3 = 2fls'm9/R, a Cartesian coordinate system (x,y,z) is 
set up so that x, y and z point eastward (the sense of spin), 
northward and radially outward, respectively. One derives 
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where /* = //(l + C A /c 2 ) is the modified j —parameter, 
(3* = [3/(1 + C\/c 2 ) is the modified (3— parameter, and 
c* 2 = (gh + C' A )/(l + C A /c 2 ) defines a characteristic speed 
c*. For free MHD tidal waves of small y— variation scales 
and with vg cx exp(iujt — ik x x — ik y y) where k x and k y are 
x— and y— component wavcnumbcrs, the cubic dispersion 
relation is to 3 - [f* 2 + c* 2 (k 2 x + k 2 y )]uj - (3*c* 2 k x = 0. Global 
versions of KARP waves are contained in eqn (5). 



3. EQUATORIALLY TRAPPED MHD TIDAL WAVES 



By a latitudinal variation of / = (3y, the equato- 
rial (3— plane becomes an effective waveguide for trapping, 
within the equatorial zone, MHD tidal waves that propa- 
gate in the </>— direction and this equatorial trapping effect 
becomes stronger in the regime of larger 4R 2 il 2 /[(gh + 
C\){\ + C 2 A /c 2 )] in eqn (5) (Longuet-Higgins 1968). We 
examine main properties of such equatorially trapped MHD 
tidal waves below. 

For equatorially trapped Alfven- Kelvin waves (Thomp- 
son 1879; Alfven 1942) around 9 = tt/2 that propagate 
only eastward at speed c*, one must directly solve eqns 
(1) — (4) with vg = to obtain 

7?, b r } = {c*,h, B} exp[-/3*2/ 2 /(2c*)]£(x - c*t), (7) 

where (/(•••) is a small arbitrary dimcnsionless waveform 
and (3* appears in the Gaussian y— profile that decays at 
large \y\ (e.g., Gill 1982). 

With / = (3y and a nonzero vg oc exp(ik x x — iut) in 
the equatorial zone around 9 = ir/2, eqn (6) becomes the 
Schrodinger equation of a harmonic oscillator. The nor- 
malized real solution takes the form of 
(3*y 2 - 
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with H n (£) being a Hermite polynomial of order n and the 
normalization constant TV = 2-"/ 2 (n!)- 1 / 2 [f3* /(ire*)} 1 / 4 . 
Variables v^, rj, and b r can be derived accordingly. The 
resulting dimensionless cubic dispersion relation for equa- 
torially trapped KARP waves becomes 

Cu 3 - (k 2 x + n + l/2)u> - k x /2 = , (9) 

where Q = uo/ {2f3*c*) 1 > 2 , k x = k x /{2(3* /c*) 1 ' 2 , and integer 
n > gives the number of nodes along y— direction. Re- 
sults of eqn (9) are displayed in Figure 1 with the straight 
dashed line for the Alfven-Kelvin mode (7) added. For 
n > 1, there are two classes of MHD tidal waves. The high- 
frequency class (upper solid lines) corresponds to equatori- 
ally trapped Alfven-Poincare waves (Poincare 1910; Alfven 
1942) that propagate either eastward or westward but with 
different phase speeds u/k x for a given \k x \. The low- 
frequency class (lower solid lines) corresponds to equato- 
rially trapped Alfven- Rossby waves (Rossby et al. 1939; 
Alfven 1942) with westward phase speeds; the group veloc- 
ity du/dk x is eastward and westward for short- and long- 
waves, respectively. The n — case (dash-dot line) of 
eqn (9) corresponds to u> = k x /2 ± (fc 2 /4 + 1/2) 1 / 2 and 
is referred to as equatorially trapped mixed Alfven- Rossby- 
Poincare waves (Poincare 1910; Rossby et al. 1939; Alfven 
1942) whose phase velocity can be to the east (plus sign) 
or west (minus sign) yet with a group velocity always to 
the east; for large k x > it behaves more like Alfven- 
Poincare waves (plus sign), while for large \k x \ with k x < 
it behaves more like Alfven- Rossby waves (minus sign). 
There are no standing modes (i.e., k x = 0) for equatori- 
ally trapped Alfven-Kelvin and Alfven- Rossby waves (Fig. 
1). In contrast, for equatorially trapped Alfven-Poincare 
waves and mixed Alfven-Rossby-Poincare waves, axisym- 
metric standing modes exist with k x = in eq. (9) and 
thus u) 2 = n + 1/2 where n = 0, 1, 2, • • •. 

4. ESTIMATES AND DIAGNOSTICS 

The preceding analysis is succinct with potential ap- 
plications to diverse stellar settings. The estimates here 
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focus on typical neutron stars. One can extend the anal- 
ysis to white dwarfs with different parameters. With 
M <~ 1.4M Q and R ~ 10 6 cm, the surface gravity of 
a neutron star is g ~ 10 14 cm s~ 2 . The ocean surface 
gravity wave speed is (gh) 1 / 2 ~ 10 9 /i 4 ^ 2 cm s _1 and the 
Alfven speed is C A ~ 10 9 .Bi 2/ o~ 1/2 cm s" 1 for HMXBs 
or Ca ~ 10 6 B 9 p 5 '' 2 cm s _1 for LMXBs, where h$, 
and B a stand for h, p, and B in units of lO^cm, 10 5 g, 
and 10"G, respectively. For an extreme magnetar of 
Ca ~ 10 12 -Bi5P 5 '' 2 cm s^ 1 , relativistic and quantum ef- 
fects are involved. For spin periods ranging from ~ 16ms 
to ~ 8.5s, the rotation speed VtR ranges from ~ 10 to 
~ 10 4 km s _1 . It is clear that MHD tidal waves are im- 
portant in the ocean dynamics of a spinning neutron star. 

For a magnetar or neutron star in a HMXB, the sit- 
uation of Ca > c, fii? and (gh) 1 / 2 may happen with a 
stronger B and a less dense ocean. The inclusion of the 
displacement current term in eqns (1), (2), (5) and (6) is 
necessary (see the definitions of /*, (3*, and c* following 
eq [6]) and the magnetic field is dynamically significant in 
MHD tidal waves besides potential diagnostic roles. 

For a neutron star in a LMXB, Ca < c, Ca < (gh) 1 / 2 , 
and Ca — unless B g ^> 1, p 5 <C 1 and h 2 -C 1. Thus, 
B affects only the rotational contribution to the speeds and 
frequencies of tidal waves significantly. Flow velocities of 
KARP waves arc powerful to push magnetic fields around 
as v 2 + v 2 > b 2 /(4irp). Meanwhile, a B of - 10 8 - 10 9 G 
would be strong enough to partially guide materials from 
an accretion disk onto the magnetic poles (White et al. 
1983), leading to intense X-ray emissions as well as spo- 
radic bursts (Bildsten & Strohmayer 1999). 

KARP waves might be driven to excitation by radiative 
thermal or nuclear processes, by polar accretion impacts, 
by magnetic thrust and dragging, or by tidal interactions 
of the accretion disk, and in turn, are capable of periodi- 
cally modulating thermal properties of the ocean and buf- 
feting magnetic fields in polar regions to affect the height 
and cross section of polar accretions. 

As an example of illustration, we now discuss possible di- 
agnostics for equatorial MHD tidal waves of a neutron star 
in a LMXB. In spherical geometry, k x = m/R for an inte- 
ger m and k x ~ 0.18m(M 33 hi) 1 / 4 (R^v?,)' 1 / 2 for a small 
C'a, where v 3 , M 33 , and Rq stand for v = CI/(2tt), M, and 
R in units of 10 3 Hz, 10 33 g, and 10 6 cm, respectively. The 
Rossby radius of deformation is D ross b y = (2(3* /c*)^ 1 / 2 

- 1.8 x 10 5 (M 33 /i 4 ) 1/4 ^ 1/2 cm. For i? 6 - 1, M 33 - 2.8 
and v 3 <~ 0.5, cquatorially trapped KARP waves influence 
a latitude range > ± 20° across the equator. 

Take a plausible geometry favorable for a detection. The 
spin and magnetic axes are misaligned with the two mag- 
netic poles at latitudes outside ±45°. The spin axis, likely 
aligned with the binary orbital axis, is taken to be perpen- 
dicular to the line of sight. During the phase of moderate 
and steady polar accretions from the surrounding disk, the 
stellar surface gives off a largely isotropic X-ray flux suf- 
ficiently away from the two magnetic polar hot spots. In 
this scenario, one may not directly detect either X-ray hot 
spots. Instead, accretion-powered KARP waves trapped in 
the equatorial belt may give rise to relatively smooth large- 
scale (L > R) periodical modulations of X-ray brightness. 
In principle, all these waves may coexist. Estimates are 



given below for /14 ~ 1, i? 6 ~ 1, and M 33 ~ 2.8. 

For m — or k x = 0, dispersion relation (9) gives stand- 
ing MHD tidal wave frequencies / w as 
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Given a detected / w , the stellar spin 
rate of v can be inferred for n — - the special case of 
mixed Alfven-Rossby-Poincare standing waves (the dash- 
dot line in Fig. 1 with k x = 0). Or, for standing Alfven- 
Poincare waves of larger n > 1, v would be smaller. 

For nonaxisymmctric waves, the frequencies of cquatori- 
ally trapped Alfven-Kelvin waves (7) (the dashed straight 
line of Fig. 1) seen in an inertial reference frame are 

\tfmv 3 



1.3 x 10 2 m(M 33 h 4 ) 1/2 R 6 2 Hz (11) 

for a small C'a, where the corotation effect is included. 
With m = 1 and a detected / w , one can infer v. For 
m = 2 and 3, / w would be precisely 2 and 3 times the 
fundamental frequency. 

For small k x - 0.18m(M 33 ft, 4 ) 1 / 4 (i? 6 ^ 3 )- 1 / 2 with m ^ 
0, the frequencies of equatorially trapped Alfven-Rossby 
waves (lower solid lines of Fig. 1) as seen in an inertial 
reference frame are 

/ w ^ 10 3 m^ 3 -1.3xl0 2 m(M 33 /i 4 ) 1/2 J R f r 2 /( 2 « + 1 ) • ( 12 ) 
With m = 1, n = 1, and a detected / w , one can infer v. 
For m = 2 and 3, / w would be approximately 2 and 3 times 
the fundamental frequency. 

If only one apparent / w were detected, it would be 
difficult to distinguish different KARP wave possibilities. 
A detection of two or more harmonics at 2 and 3 or 
more times the fundamental frequency may offer valu- 
able hints. For nonaxisymmetric equatorially trapped 
Alfven-Poincare waves and mixed Alfven-Rossby-Poincare 
waves traveling around the globe, relation (10) must be 
modified to include the corotation effect mv and other 
m— dependence (Fig. 1); the overall m— dependence of 
these frequencies is not linear. Similarly, the more accu- 
rate version of relation (12) for Alfven-Rossby waves does 
not contain a strictly linear to— dependence either (Fig. 1). 
By detecting more harmonics, one can constrain the type 
of KARP waves involved according to the dispersion rela- 
tions (Fig. 1). As f w given by eqns (10) — (12) vary with v, 
the expected stellar v spin-up in a LMXB (or unexpected 
variations, such as accreting X-ray pulsar Cen X-3) would 
lead to a / w spin-up (or corresponding variations) . 

There have been recent theoretical development on stel- 
lar r— mode (Papaloizou & Pringle 1978; Saio 1982) insta- 
bilities caused by gravitational radiations (Andersson et 
al. 1999) to limit spin rates of hot young neutron stars 
at birth. The prospect that fast spinning neutron stars 
in LMXBs might emit gravitational waves (Bildsten 1999) 
observable by ground-based detectors currently under con- 
struction (e.g., the enhanced LIGO) offers another poten- 
tial test for the KARP wave scenario advanced here by 
independent determinations of neutron star spin rates. 

In the case that the stellar surface of a radio pulsar is 
covered with an ocean, various MHD tidal waves may mod- 
ulate or interfere with radio pulse emissions from the mag- 
netic polar regions by (quasi-)periodically varying polar 
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magnetic field strengths or changing magnetospheric con- 
ditions (Lou 2001a). This may lead to observable effects 
because magnetic and spin axes misalign for radio pulsars 
to shine. With B > 10 12 G for a pulsar, C A and {gh) 1 / 2 
can be comparable and be both > ViR. For equatori- 
ally trapped Alfven-Poincare waves, Alfven-Kelvin waves, 
and the high-frequency branch of mixed Alfvcn-Rossby- 
Poincare waves, modulation timcscales would be shorter 
than or comparable to radio pulse periods. For equatori- 
ally trapped Alfvcn-Rossby waves and the low-frequency 
branch of mix Alfven-Rossby-Poincare waves, modulation 
timescales would be longer than radio pulse periods. These 
are intriguing possibilities to be explored further. 
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